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Abstract. We present detailed photometric investigations of the recurrent nova CI Aql. New data obtained after 
the 2000 outburst are used to derive a 3D geometrical model of the system. The resulting light curves clearly 
indicate the existence of an asymmetric spray around the accretion disk, as claimed in the past e.g. for the 
super-soft X-ray source CAL87 in the LMC. The simulated light curves give us the mass transfer rates vary- 
ing from M ^ 2.5 x 10"* Mq yr"^ in 1991-1996 to 5.5 x 10"** < M < 1.5 x 10"'' Mgyr"^ in 2001/2002. The 
distance and the interstellar foreground extinction resulting from the model are 1.55 kpc and Eb-v ~ 0™98 re- 
spectively. During fast photometry sequences in 2002 short timescale variations {tp « 13 minutes) of the mass 
loss are found. Moreover a change in the orbital period of the system is detectable and results in a mass loss of 
2.2 X 10"*^ < AM < 5.7 x 10"® M© during the nova explosion. 
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1. Introduction 

CI Aql is one of the 9 to 10 members of the class of recur- 
rent novae: U Sco, V394 CrA, RS Oph, T CrB, V745 Sco, 
¥3890 Sgr, T P yx, C I Aql and IM N or (W ebbink et a l. 
19871; Sekiguchi [1995^ Schmeja et al. |2000t Liller |2002D . 



Webbink et al. (1987) also mention V1017 Sgr as possi- 
ble class members although its status is still not clear. 
The first known outburst of CI Aql was discove red o n 
Heidelber g plat es recorded in June 1917 (Reinmuth |l 925|) . 
Williams ( ^000|) completed the light curve by using records 
on Har vard College Observatory patrol plates. Schaefer 
(2001a) foimd another outburst in 1941 again on Harvard 
plates. Schaefer argues, it might be a recurrent nova with 
a timescale of 20 years and the 1960 and 1980 outburst 
were missed. As the timescales of other recurrent novae 
often change and as there are no observations available, 
we assume for our calculations a quiescence phase of 60 
years before the 2000 event. Anyway this does not affect 
the results of the model presented here but our resulting 
pre-outburst accretion rate indicates a long recurrence 
timescale. CI Aql was found to be an eclipsing binary 
system with a period of 0^618355(9) by Mennickent & 
Honeycutt (1995). It is, to our knowledge, the only eclips- 
ing system investigated photometrically in such detail be- 
fore and after an outburst. Following the classification 



of Sekiguchi (199£), CI Aql is of U Sco subclass (shghtly 



evolved main sequence star and accreting white dwarf). 
For a detailed disc ussion of the outburst data we refer to 
Matsumoto et al. (|200l|) and Kiss et al. ( ^00 1[) . 
In this paper we deduce a detailed model of the sys- 
tem basing on our optical photometry of 2001 and 2002. 
Further we follow the final decline to quiescence and find 
quasi periodic short timescale variations in thi s syst em. 
Together with pre-outburst data of Honeycutt (2001) we 
finally determine a period change. 



2. The Data 

The new data were obtained with the Innsbruck 60cm 



telescope (Kimeswenger 2001) and a direct imaging CCD 
device in the period from June 21, 2001 to July 9, 2002. 
In 2001 a CompuScope Kodak 0400 CCD (4f6 x 3.'1 field 
of view) was attached, in 2002 it was an AP7p SITe 502e 
(8.'36 X 8f36). 781 images were taken in 34 nights with 
V, R and Ic filters. The exposure times varied with filter, 
brightness, weather conditions and camera. Flatfield and 
bias subtraction were carried out in standard manner with 
the help of MIDAS routines. The source extraction was 
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performed using SExtractor V2 (Bertin & Arnouts 1996). 
The rms of the comparison standards in the field was 
< 0'.°08 with the Kodak chip and < 0'.°05 with the SITe 
CCD typically. The light curve was obtained by means of 
differential photometry of up to 40 stars within about 4.'0 
from the target in this very crowded field. For the abso- 
lute calibration CCD standards in the field by Henden & 
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Honeycutt ( |1995D and Henden ( 2001 ) were applied. We 
found on that occasion that the coordinates of the whole 
set of the Henden & Honeycutt standards are shifted 1.'73 
east and 1.'20 south and thus the finding chart overlay 
(SIMBAD/ALADIN) is wrong. As w e use Johnson R and 
the CCD standards of Henden (2001) are taken in i?c we 
have to assume color terms in the absolute calibration. 
The attempt to derive these terms failed within the accu- 
racy of the photometric data. We thus did not use the R 
band for the final results in extinction and distance (see 
section ||). 

The pre-outburst data ranging from June 4, 1991 to 
Septem ber 29, 1996 were provided to us by Honeycutt 
(2001). The data consist of two sets (changing June 1995) 
with different zero points. The first set (the one used for 
Mennickent & Honeycutt 1995| ) was shifted according to 
the information given by Honeycutt to overlay. This cor- 
responds well to the calibration of Skody & Howell ( 19921 ). 
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3. The Photometric Behavior 

First of all we derived the 'remnant' of the outburst in our 
2001 data. The object clearly had no t return ed completely 
to quiescence as stated by Schaefer ( 2001b|) for August 4, 
2001. To verify the real level, the new period (see section ^ 
was taken and the data points outside the primary eclipse 
were used to derive the final decline of the nova outburst. 
We see (Fig. |^) that this decline was finished somewhen 
in February/March 2002. 
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Fig. 1. The final decline phase of the recent nova outburst 
(upper: /c; lower: V). Data around primary minima were 
removed to decrease effects due to the orbital modulation. 
First order polynomials to the 2001 data and second order 
polynomials to the whole data set are shown. 

In contrast to this slowly continuous decline, the 



plateau phase in the model of Hachisu & Kato (2001) 
should stop abruptly when the super-soft X-ray source 
phase ends. The position for this end of the plateau 
phase depends strongly on the hydrogen content X. The 
overall evolution is shown in Fig. |[ We find that the 
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Fig. 2. The complete decline phase of the 2000 outburst 
(open symbols: CCD data from VSNET; closed symbols: 
our photometry) . Due to the orbital modulation the upper 
boundaries have to be used. The ticks mark the predicted 
dates for the abrupt decline after the plateau for different 
hydrogen content X (Hachisu & Kato 2001). 



decline started between JD 2452050 and JD 2452100. 
This results, including th e late decline model calcula- 
tions of Hachisu & Kato (2001), in a hydrogen content 
of 0.45 < X < 0.50. This is in contradiction to Hachisu & 



Kato (2002), who obtained a hydrogen poor model with 
X = 0.35 on the basis of only two data points from the 



lAU Circulars of Schaefer (|2001b|, |2001c|) . Our moderately 
hydrogen rich result is supported by the infrared spec- 
troscopy of Lynch et al. ( |2002 ) and by the first decline 
model (Fig. 4 in Hachisu & Kato |200lD . 
The linear decline (Fig. |l|) was subtracted to derive the 
light curve for 2001 which is shown together with those 
from the pre-outburst and from 2002 in Fig. 3. The peri- 
ods and epochs used are described in section 7. 

4. The 3D Model 

The double star system was modelled with a slightly 
evolved main sequence secondary (SE) filling the Roche 
lobe, a nearly point source like, spherically symmetric 
white dwarf (WD), a rotationally symmetric inner accre- 
tion disk, a thin accretion stream from the Lagrange point 
towards the disk and a spray of bounced material first in- 
troduced by Schandl et al. (|1997| ). 

Th e WD mass of 1.2 Mq was derived by Hachisu & 
Kato ( 2001 , Fig. 4) by means of the thick wind model of 
the early decline of the 1917 and the 2000 outburst. The 
mass of the SE can not be determined that directly. For 
a range of SE masses we used the evolutionary tracks and 
colors for solar abundance stars of Girardi et al. (200C ) and 
our 2002 photometry at minimum, when the SE dominates 
the emission. The point in the evolutionary track was cho- 
sen to fill the Roche lobe. This results in the extinction 
free (t^ — /c)o color. Together with our measured {V — Ic) 
in 2002 and each SE mass we are able to derive a corre- 
sponding range for Eb-v- The results are summarized in 
Tab. |l|. 
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Fig. 3. The V light curves for three different epochs: 1991- 
1996 (upper), 2001 (middle) and 2002 (lower panel). The 
light curve for 2001 was corrected according to the lin- 
ear final decline shown above and normalized relative to 
October 1, 2001. The scatter is dominated by the short 
timescale variations (see section ^) and not by the errors 
in the data. 

Table 1. Extinction derived from {V — Ic)minimum = 1™5 
and the evolutionary tracks at Roche lobe filling position. 
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The geometric model was realized with the help of 
t he ''M ATLAB language of technical computing 5.3R1V' 
( |1999D . This allowed us to easily implement and modify 
various asymmetric surfaces, system inclinations and rota- 
tion periods. All surfaces are assumed to radiate like black 
bodies. The ray tracing for the irradiation within the sys- 
tem was solved explicitly. The ray tracings during system 
rotation for the resulting light curves were obtained by 
the internal renderer. The model consists of the five com- 




Fig. 4. The basic geometric model with the Roche lobe 
filling SE, the symmetric accretion disk, the accretion 
stream and the spray of bounced material at the light 
curve phase of 0.33. The grid mesh density of the sur- 
faces was reduced here for clarity by a factor of three with 
respect to the one used in the real calculations. 



ponents described in detail later. The resulting scheme is 
shown in Fig. ^. Each surface element was characterized 
by its temperature. To obtain the flux in each band, the 
contribution of each element was calculated by folding the 
black body with the filter curves {B, V, Rc and Iq- Bessel 
AUer et al. 



1990; Johnson R: AUer et al. 1982). For the zero point a 
spherical body simulating the sun was used in the same 
way in the renderer and calibrated by the values of Bessel 
et al. ( |1998| ). 



The White Dwarf 

We assume a spherical WD photosphere. The luminosity 
of the WD is given by 



According to van den Heuvel et al. ( 1992 ) a 1.4 M© 
ZAMS star as SE can already achieve a mass trans- 
fer, which leads to a steady hydrogen shell burning on 
the WD. A recurrent nova is not likely anymore below 
this mass limit. Thus a smaller SE does not have to be 
taken into account. On the other hand the VI^JK^ pho- 



tometry of Schmeja et al. (200C) shortly after the out- 
burst limits the interstellar extinction. The {V — Ks)q 
as well as the (/c — J)o were calculated for different in- 
terstellar extinctions. They get to a negative domain at 
Eb-v > 1™05. This is thus an upper limit for the extinc- 
tion, assuming a dominant photosphere. Thus we conclude 
1.4 < MsE < 1.7 Mq. 



1 GA/wpM 

2 Rwu 



(1) 



where Xhsb is the luminosity obtained from the remnant 
of the hydrogen shell burning in 2001 (Hachisu & Kato 
^001 ) and the second term originates from the accretion 
(Shaviv & Starrfield 1987). According to the model of 
Hachisu & Kato and the photometric behavior shown be- 
fore, there is a remnant of the thermonuclear processes 
after the outburst in the 2001 data. The simulations give 
us Lhsb ~ 200 L0. This results in a net hydrogen mass 
burning A/hc ~ 10~^ Mq during the plateau phase from 
end 2000 to August 2001. This is just below Vio of the 
mass-increase rate given by Hachisu & Kato for the total 
outburst. 
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The direct light contribution of the WD is negligibly 
small in our visual bands. It only contributes by the irra- 
diation on the other components. 

The Secondary 

The shape of the secondary is calculated numerically. It 
starts in the inner critical Lagrange point which is fixed 

by 



(2) 



where G is the gravitational constant and Twd, ^se and 
the vectors to the WD, the SE and the center of 
gravity (CG) respectively, is the unity vector parallel 
to J'cG ~ ''sE- An index stands for values at the Lagrange 
point. The surface of the SE then follows the equipoten- 
tial (p — (pQ including the rotation of the complete binary 
flattening the system 




(3) 



The calculated volume corresponds within one percent to 
the Roche lobe volume of Eggleton ( 1983 ). 
The star is irradiated by the WD. This increases the undis- 
turbed temperature T^,^: 



47^7-2 



COS ( 



(4) 



where r is the distance of the surface element from the 
WD and the angle of incidence. An efficiency rj^, — 0.7 
was used. The value results from fitting details of the light 
curve - mainly the gradient at phases 0.08 to 0.18. This 
somewhat higher value than that assumed in previous in- 
vestigations may originate from the compact MS nature 
of the SE preventing the star from using the irradiated 
energy to change its internal structure. For the model this 
energy is distributed along the stellar surfac e (Fig . ||) in 
the same way as described in Schandl et al. (1997). 




meridian declination [ 



Fig. 5. The effects of the redistribution of the temperature 
along the meridian from the inner critical Lagrange point 
(0°) to a pole (90°) and further to the backside of the SE 
(180°): thick line - redistributed temperature profile; thin 
line - original irradiated temperature T^ i^r- 



The Disk 

A rotational symmetric accretion disk was used. 



Elongated disks as they are proposed by Paczyiiski (1977) 
or result from projecting the 3D disk models of Hirose et 
al. (1991) were tested too but showed neghgible effects. 
The horizontal size i?D = 0.75i?^n of the disk is deter- 
mined by the stability due to tidal forces near the Roche 



lobe (Paczynsky 1977, Papaloizou & Pringle 1977). The 
width of the primary minimum independently leads us 
to the same result. The vertical height in the inner part 
(r < Ri) was obtained from the fit to the hydrodynamic 
models (Meyer & Meyer-Hofmeister 1982) by Schandl et 



a.1. (1997 ). The calculations of Meyer & Meyer-Hofmeister 
( 1982| ) show a rapid turnover of the vertical behavior. The 
lower gradient causes for those regions no irradiation by 
the WD anymore. Therefore we used 

2.527 ri 093 A/-0.38 ^.;^o.i7 . ^ < 
z = ^ (5) 
const. : Ri < r < i?D 

where z, r, Mwd and M are in solar units. i?i was es- 



timated from Fig. 1 in Meyer & Meyer-Hofmeister (1982) 
for M„D = 1.2 Mq to be 



Ri 0.365 log(M) + 2.69 



(6) 



This description of the disk differs significantly from the 
parabolic law used by Hachisu & Kato ( |200l| ). They ar- 
gue, referring to Schandl et al. (1997), that they simulate 
the "flare-up" of the rim. But they increase strongly the 
reradiation of the WD flux on every part of the disk. This 
effect gets drastically at the outer disk portions. Moreover 
the high "flare-up" built by the spray, as shown below, 
does not cover the whole rim but only about a quarter of 



it. Meyer & Meyer-Hofmeister (1982) also tested the ef- 



fects of the irradiation on the vertical structure and found 
only marginal effects for CVs. 

The disk temperature was adopted from the frictional 
heating and the surface irradiation as in Schandl et al. 



( I1997D by 



4 3 GAfvs 

= Z ■■ 

air r 



M 



47^7-2 



cosf 



(7) 



For Td < TsT the stream temperature Tst (see descrip- 
tion of the accretion stream later in this section) of the 
infalling material was assumed at the disk too. The effi- 
ciency parameter for the irradiation was assumed to be 
rjD = 0.5, like in Schandl ct al. (1997) and Hachisu & 
Kato ( ^001 ). Our tests with the models show nearly no 
noticeable effects on the composite light curve for slightly 
different rjjj . The black body approximation fits very well 
for the disk in the visible spectral range as shown by recent 
full radiative transfer models (AcDc) producing artificial 



spectra (Nagel et al. 2002) 



The Spray 

The spray was first introduced by Schandl et al. (1997) 



as convolute of individual blobs bounced after reaching 
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Fig. 6. The model for the different phases, 
distribution. 



The grey scale represents the V flux intensity due to the temperature 



t he in cident point of the accretion str eam. S chand l et al. 
(|1997D and Meyer-Hofmeister et al. ( |l997t |l99S|) found 
that the vertical extension dominates the effect on the 
light curve for the high inclination system here. Although 
they were able to calculate individual trajectories of undis- 
turbed blobs the final, optically thick surface was mod- 
elled from the effects in the light curve (see Fig. We 
obtained the shape of the spray in the same way for our 
target (Fig. 0). The massiveness of the spray varied with 
the mass transfer rates. 



The temperature of the irradiated regions of the spray 
were deduced according to Eqn. ^ with the stream temper- 
ature as undisturbed value. Since the material consists 
of a clumpy medium the overall efficiency for heating due 
to irradiation is not comparable to a stellar photosphere. 
At short wavelength it i s assu med that the spray is even 
half transparent (Meyer 2002 ) causing the observed super 
soft X-ray effects. We therefore found in the simulations 
the plausible value of r]s = 0.2. 



The Accretion Stream 

The end points of the accretion stream are defined on the 
one hand by the Lagrange point and on the other hand by 
the impact point. The latter is calculated from an undis- 
turbed gravitational particle trajectory towards the disk 
and characterized by the angle (j)^ . As the bend of the tra- 
jectory is small (Fig. 0), a straight line was assumed. The 
cross-section of the stream is determined like in Meyer 
& Meyer-Hofmeister ( |l983 ) but slightly flattened. Recent 
hydrodynamic numerical simulations by Oka et al. (2002) 
show very similar deflection angles and tube geometries at 
r w 1. 



The temperature Tst of the material is defined by the 
irradiated and distributed temperature of the SE at the 



stream source. The contribution of the stream to the light 
curve is rather small. 



The Composition of the Modell 

The individual components are dominating different parts 
of the light curve: 

— The width of the central part of the primary minimum 
gives the disk size and limits the system inclination i. 

— The depth of the primary and the secondary minimum 
strongly depend on the system inclination. Neither an 
inclination of 70° nor of 72° result in a usable fit, as 
in the first case the modelled primary minima are not 
deep enough for reasonable M and adjusted secondary 
minima, otherwise those of 2002 are too deep. 

— The difference of the plateau levels before and after the 
secondary minimum gives the parameters of the sec- 
ondary - namely the reduced irradiation by the shadow 
of the spray - and the irradiation efficiency of the spray. 

— The depression of the late light curve gives length and 
height of the spray. 

— The gradient of the late part of the primary minimum 
(phase: 0.08 - 0.18) gives the irradiation parameter r/sE 
of the SE 



The parameters Mse, Tse, h M, Rd, 77,, rju, rjs and 
the spray geometry were varied independently in the be- 
ginning. The use of different photometric bands - namely 
V, R and Ic - gives us further restrictions to the over- 
lapping part of the possible individual solutions. Finally 
MsE, Tse, i, V*, Vd and 775 are limiting the variation 
ranges between the years. The pre-outburst parameters 
are missing other bands than V and thus are somewhat 
more flexible. The resulting simulated light curves are 
shown together with the data in Fig. ||. In Fig. ||the indi- 
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Fig. 7. The structure of the spray as top view (upper 
panel, including the circular disk) and its vertical struc- 
ture projection (lower panel). The shape is symmetric with 
respect to the equatorial plane. The coordinates originate 
in the center of gravity (GC). For the stream the real 
trajectory is shown here. For the calculation a straight 
approximation was used. 

Table 2. The parameters of the best fit model. 



system: 

secondary: 

white dwarf: 

disk: 

stream: 

1991-1996 

2001 

2002 



i = 71° 

MsE = 1.5 M© 
MwD = 1.2 M© 
Rd = 0.75 -RwD 
<^s =14? 2 

M = 2.5 X 10"** Mq 
M = 1.5 X 10"^ 
M = 5.5 X 10"* 



TsE = 7000 K 
i?wD = 0.0072 Rq 



Ri =0.54i?^'^> 



Ri = 0.75 
Ri = 0.72 



vidual contributions of the components to the total light 
curve are plotted. 



5. Distance and Interstellar Extinction 

The model light curves were derived independently for 
each band. When shifting those curves to the data we ob- 
tain for each band an optimized solution as function of the 
distance D and the foreground extinction. Having more 
than two bands allows us to derive Eb-v (Fig. |l^), using 




phase 



phase 



Fig. 9. The contribution of the individual components to 
the total hght curve in 2001 for V (left) and Ic (right). 
The asymmetry of the spray causes a shift of the position 
in the primary eclipse. This is of special interest if combin- 
ing data from different epochs (respectively mass transfer 
rates) for determination of the period. 



a standard model for the interstellar extinction (Mathis 
et al. |1977]). 



Kiss et al. ( pOOlD derive values of 0'?"83±0'.°20, 
1^"08±0'.°20 and 0'."66±0?"30 for the diffuse interstellar 
bands (DIB) at 584.9 nm, 661.3 nm and Call 393. 4 nm 
respectively. Usi ng the photometries of Hanzl (200C) and 
Jesacher et al. (200C) on days +2 and -f6 after maxi- 
mum and the general tendency for novae around maxi- 
mum of (B-V)o a"23±Q^06 (Kiss et al. |200lD we get 
0?^35 <Eb_v < 0™75. But, as Kiss et al. point out al- 
ready, those colors may be affected strongly by emission 



features. Hachisu & Kato (2001) derive from their model 

-V = 0'.°86. In their re- 



of the pre-out burst light curve En- 
vision (Hachisu & Kato 2002D , based on the early final B 
band decline points of Schaefer (2001b, 2001c) and thus 
the low hydrogen content (see section |^), they move it 
to Eb-v = 1"^00. While the direct spectroscopic methods 
give similar results, the photometry suffers from the fact, 
that an average color for novae is assumed. The compar- 
ison with the models using the B band may also suffer 
from the fact, that this band is mostly affected by the non 
grey opacities. A black body, as used by the models, may 
not work properly. We mostly rely on the direct spectro- 
scopic methods and the models here using the red bands 
and thus use Eb-v = 0"98 ± 0'"1. 

The total luminosity also results from the simulations. 
This gives us a distance of 1.52 kpc to 1.58 kpc. Hachisu 



& Kato (2001) get 1.6 kpc from the maximum outburst 
brightness and using Eb-v = 0™85. Whereas they derive 



1.1 kpc in their revision (Hachisu & Kato 2002) 
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Fig. 10. The distance as function of the interstellar ex- 
tinction for different bands individually. The solutions 
are derived by the intersection of the pairs of each year. 
The thick lines are the solution for 2001 resulting in 
Eb-v = 1^025. The thin lines represent the 2002 data 
giving Eb-v = 0'."935. The dashed line represents the 
pre-outburst data. Due to the lacking second photometric 
band, no solution is possible there. Nevertheless the ten- 
dency gives us some information on the accuracy of the 
solutions above. 



expelled due to variations and instabilities of the mass 
transfer flow. This causes consequently a temporarily 
extended spray, detectable as small "outburst". Meyer- 
Hofmeister et al. use f ^K: where is the Keplerian 
orbital period at the edge of the accretion disk. They 
assume that the blobs oscillate free around the mean 
circular orbit and thus give the timescale for the "out- 
bursts". We, in contrast, use the angular dimensions 
of the spray. The spray will orbit with about the same 
rotational velocity as the accretion disk. The geometric 
models result in a spray, where the expelled material re- 
turns after about one quarter of an orbit to the accretion 
disk. This results in i or in our case 7™" (0.008 P). 
As the blobs of the spray reach out to up to 1.4 times the 
accretion disk radius the timescale can be up to 12.5™'" 
(0.014 P). The finite duration of the accretion outbursts 
extend those times. Thus the calculated timescales are 
the lower boundaries. In fact we find these variations 
during the minima, when the extended spray geometry is 
well visible and at the same time the luminosity contrast 
is preferring the effect in the photometry (Fig. O). 



6. Fast Mass Transfer Variations 

The bandwidth of the photometric variations, as shown in 
Fig. ^, gives us information on short timescale variations 
of the mass transfer. Meyer-Hofmeister et al. (1998) found 
already in RX J0019.8-)-2156 variations in the timescale 
of 1^2. To monitor the details of the mass transfer we 
carried out fast Iq photometry in the nights of June 26, 
2002 and July 4, 2002. Both nights were characterized 
by a primary minimum. The minima were overlayed and 
the lower boundary was used to define the undisturbed 
minimum (Fig. pi] ). 



14.4 



14.6 



14.8 



15.0 




-0.1 



-0.05 





phase 



0.05 



0.1 



Fig. 11. The Ic minima of June 26, 2002 (circles) and 
July 4, 2002 (squares). The hne corresponds to the lower 
boundary limit fit. 



After the substraction 
ima, the residual show 

0.011 < < 0.018 P (10' 

Hofmeister et al. (|1998D 



of the undisturbed min- 
the typical timescale of 
< < 16""'"). Meyer- 
assume that some of the 
individual blobs forming the spray are more violently 
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-0.02 



0.00 0.02 
phase 



0.04 
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Fig. 12. The residuals of the Ic minima of June 26, 2002 
(upper) and July 4, 2002 (lower). The lines are two pe- 
riod floating means. The timescale of the small outbursts 
corresponds to 10""'" < ip < 16™'". 



7. The Period - The Outburst Mass 

Our data pointed out a shift of the period. The periods 
were calculated by using the PDM method (Stellingwerf 
1978). The cores of the 9 minima (Fig. |l^) were fitted 
by polynomials to define the exact minima. To obtain ho- 
mogeneity, the data by Honeycutt ( |2001| ) were used to 
recalculate the pre-outburst parameters. To demonstrate 
the necessity of this recalculation Fig. [l^ shows the PDM 
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minima, containing only the first data set by Mennickent 
& Honeycutt with the default sampling parameters of 
Stellingwerf (197S). The absolutely lowest noise peak leads 
to the original result. Recalculation with optimized pa- 
rameters for PDM shows different results. We thus pro- 
pose the light curve elements of 

2448412.147(46) 0':'6183609(9) x 

for the pre-outburst phase. For the post-outburst phase 
we get 

2452081.5022(46) + 0^6183634(3) x E. 



Note that the epoch for the minima will slightly change 
with decreasing accretion rate (Fig. ||). 




0.6183 



period [d] 



0.6184 



Fig. 13. The PDM 6 cores of the V data around the min- 
imum. The thick lines give the results from the data (up- 
per: Honneycutt 2001 ; lower: our data). The grey lines 
are polynomial fits of the minima cores to derive the ex- 
act positions. The arrows mark the derived minima and 
t he or iginal period derived by Mennickent & Honeycutt 
(119951) . 



The period above is supported independently by the /c 
data giving us a period of 0'?6183627 (= 0.06 seconds less). 
Using the change of the period P 



2„3 



P' = 



GM 



{P + APf = 



47r^a' 



2-3 



G (M - AM) 



(8) 



where a and a are the separations before and after the 
outburst, M the total mass of the system, AP the period 
change and AM the outburst mass. The conservation of 
angular momentum gives 




0.6182 



0.6183 0.6184 
period 



0.6185 



Fig. 14. The PDM diagram of the original data set of 
Mennickent & Honeycutt (1995) covering 1991 to April 
1995 only by using the default parameters for the binning 
by Stellingwerf (1978). The original result is indicated at 
the position of the lowest noise peak. 



a 



[ MwD 



( MSE , 



p 



AM 



MsE 



(9) 



a2 



P + AP 



AM 



A/wD-AAf 



1 



The parameter G [0, 1] refers to the efficiency of the an- 
gular momentum transfer to the expelled mass. The two 
extreme are that the expelled mass carries away its origi- 
nal angular momentum (^ = 1) and the opposite case that 
all angular momentum is transferred and thus kept in the 
system due to viscosity (f — 0). We use these two cases 
for the calculated range of mass loss given below. The re- 
sulting expelled mass does not depend on the period P 
within the tested range (100 times the formal error). Also 
the mass of the secondary (testing 1.2 <Mse < 2.5 Mq) 
does not affect the result. The mass of the WD scales 
direct proportional to AM. As it is known to a 10-15% 
level, it does not influence AM significantly. The main un- 
certainties originate from the errors of AP and from the 
choice of ^. Within a factor of 10 around the determined 
value we got AP oc AM. This gives us, using the extreme 
values for ^ and three times the formal error of AP 

2.2 X 10"^ MfT, < AM < 5.7 x 10"^ M^. 



Hachisu & Kato ( |2001| ) obtained 4.7 x 10"^ Mq from 
the thick wind model of the outburst. Assuming the 
last outburst 60 year ago, our result corresponds 
to 3.6 X 10'^ Mo/yr < M < 9.5 x 10"^ Mg/yr. This 
corresponds well to the accretion rates of the models dur- 
ing the quiescence phase. 
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8. Conclusions 

The models presented allow us to determine the physical 
parameters of the system within a very small range. This 
is mainly achieved by the combination of data from dif- 
ferent bands and different epochs what is the main ad- 



vantage to the work of Hachisu & Kato (2001 



The model of the recurrent nova two years after the out- 
burst, when it reached a mean quiescent photometric state 
still differs significantly from the one of the pre-outburst 
phase. Even speculative, one may assume that the energy 
transferred during the outburst, when the secondary was 
completely enclosed by the WD shell, caused a signifi- 
cant extension over the Roche lobe. This increased the 
mass transfer to 5.5 x 10"^ < Af < 1.5 x 10""^ Moyr"^ 
in 2001/2002. Thus if the mass transfer rate derived from 
the pre-outburst phase is used as average, it does not 
reach up to the accreted mass during the 60 years of quies- 
cence. The mass transfer rate of M ~ 2.5 x 10~^ Mq yr~^ 
derived for 1991-1996 differs from that given by Hachisu 
& Kato by a factor of four due to the different model 
of the accreting region. The period shift determined in 
this paper gives us a good estimate of the expelled mass 
2.2 X 10"^ < AM < 5.7 x 10"^ Mq. Within the errors 
(the uncertainty due to the angular momentum carried 
by the ejected material and the uncertainty of the aver- 
age mass transfer throughout the 60 years of quiescence) 
the determined expelled mass corresponds to the accre- 
tion. Thus the evolution towards the critical mass of the 
WD runs rather slow at this evolutionary stage. Assuming 
a net mass increase of < 10~^ Mq per outburst leads to 
a few 10^ years to reach the critical mass. On the other 
hand the evolutionary tracks of the secondary show it at 
a position with extremely rapid increase of its diameter 
(about 5 10"^^° Rq yr~^). Moreover the decrease of the or- 
bit is about 10"'' R© per outburst. Thus the mass transfer 
should increase and may even evolve towards steady hy- 
drogen burning (van den Heuvel et al. 1992 ). The derived 
distance and interstellar extinction gives us a somewhat 
higher luminosity than in the outburst models of Hachisu 
& Kat o (^001 ). This affects also the considerations in Oka 
et al. ( 2002| ). Because of our smaller inclination we also 
see some parts of the disk during the primary minimum 
and thus it is not necessary to increase the temperature 
of the SE in 2001 to obtain a higher flux. 
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